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We develop a finite temperature perturbation theory (beyond the mean field) for a Bose- 
condensed gas and calculate temperature-dependent damping rates and energy shifts for Bogolyubov 
excitations of any energy. The theory is generalized for the case of excitations in a spatially inho- 
mogeneous (trapped) Bose-condensed gas, where we emphasize the principal importance of inho- 
mogeneity of the condensate density profile and develop the method of calculating the self-energy 
functions. The use of the theory is demonstrated by calculating the damping rates and energy shifts 
of low-energy quasiclassical excitations, i.e. the quasiclassical excitations with energies much smaller 
than the mean field interaction between particles. In this case the boundary region of the condensate 
plays a crucial role, and the result for the damping rates and energy shifts is completely different 
from that in spatially homogeneous gases. We also analyze the frequency shifts and damping of 
sound waves in cylindrical Bose condensates and discuss the role of damping in the recent MIT 
experiment on the sound propagation. 



(N 

> . I. INTRODUCTION 

in 



Recent developments in the physics of ultra-cold gases have led to the discovery of Bose-Einstein condensation 
(BEG) in trapped clouds of alkah atoms [|l]-§ and stimulated a tremendous boost in theoretical studies of weakly 
f*^ interacting Bose gases. As in previous years, these studies rely on the binary approximation for the interparticle 
00 , interaction. The latter is characterized by the 2-body scattering length a, which assumes the presence of a small 
On . gaseous parameter no? ^ 1 (n is the gas density). Especially intensive are the attempts to reach beyond the ordinary 
mean field approach and to develop a theory which can properly describe the behavior of finite temperature elementary 
excitations of a trapped Bose-condensed gas and in particular, explain the JILA Q and MIT Q experiments on energy 
shifts and damping rates of the excitations. 

The commonly used mean field theory (for a > 0) is based on the Bogolyubov quasiparticle approach developed 
originally for a spatially homogeneous Bose-condensed gas at T ^ Q and employed by Lee and Yang Q (see also 
Q at finite temperatures. The generalization of the Bogolyubov method for spatially inhomogeneous systems has 

O . been described by De Gennes . In the case of a Bose-condensed gas it should be completed b y th e equation for the 
I wavefunction of the spatially inhomogeneous condensate, derived by Pitaevskii and Gross [ll| . 
• I— I : For spatially homogeneous gases the theory beyond the mean field approach was also developed. Beliaev p^ ] 
constructed the zero-temperature diagram technique which allows one to find corrections to the energies of Bogolyubov 
^ , excitations, proportional to (noa'^)^/^, where uq is the condensate density. The corrections are provided by the 
■ - - • interaction between the excitations (in particular, through the condensate) and contain both real (energy shift) and 
imaginary (damping rate) parts. At T = the latter originates from spontaneous decay of a given excitation (y) to 
two other excitations (7 and 7'), with smaller energies and momenta: — > 7 -f 7'. A universal expression for the 
chemical potential in terms of the self-energy functions has been found by Pines and Hugenholtz It should be 
emphasized that the corrections proportional to ni^a^ already depend on the contribution of 3-body interactions and, 
hence, can not be obtained within the binary approximation. 

The Beliaev approach was employed by Popov at finite temperatures. In this case the corrections to the 
Beliaev self-energies contain infra-red singularities, i.e. they tend to infinity for momenta p — > 0. This prompted 
Popov to make a renormalization of the theory, which links the microscopic approach with phenomenological Landau 
hydrodynamics |l5|| . The Popov theory eliminates the infra-red singularities and allows one to describe the behavior 
of low-energy excitations (phonons) at temperatures much smaller than the mean field interparticle interaction uqU 
(U = Anh^a/m, with m being the atom mass). The damping of phonons in this temperature range is determined by the 
Beliaev damping processes and has also been calculated by Hohenberg and Martin ]T6t . A simplified approach within 
the dielectric formalism was used by Szepfalusy and Kondor for calculating the damping rates of excitations in the 
phonon branch of the spectrum. They found that at temperatures T 3> tiqU the damping rate of a given excitation (v) 
originates from the scattering of thermal excitations (7 and 7') on the excitation v through the processes 1/-I-7 7'. 
Since the characteristic energies of the thermal excitations 7, 7' turn out to be much larger than the energy of the 
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excitation this damping channel can be treated as Landau damping. It should be noted that the damping rates 
can be simply found by considering the interaction between the excitations as a small perturbation and using Fermi's 
golden rule. This allows one to properly take into account the Bogolyubov nature of the thermal excitations. The 
damping rates of phonons in a spatially homogeneous Bose-condensed gas, in particular for the Szepfalusy-Kondor 
mechanism, have been calculated in the recent contributions p^-[2^] . 

In order to reach beyond the mean field theory at T ^ noU one should further develop the Popov approach. One can 
also proceed along the lines of the Beliaev theory, since any physical quantity should be determined by combinations 
of the Beliaev self-energies, which do not contain the infrared singularities. In this paper we choose the latter way 
and construct the perturbation theory for a Bose-condensed gas, which allows us to find the next to leading order 
terms (the terms proportional to {nQa^Y^'^) in the energy spectrum of the elementary excitations. As in |l7|-p2|, we 
consider the excitations in the so-called collisionless regime, where their De Broglie wavelength is much smaller than 
the mean free path of the thermal excitations. 

We start with the case of a spatially homogeneous Bose-condensed gas and find temperature-dependent energy 
shifts and damping rates for Bogolyubov excitations of any energy. At temperatures T ^ noU the small parameter 
of the theory proves to be 

^Ka^)^/^«1, (1) 

in contrast to ngo'^ ^ 1 for T ~ 0. The appearance of the extra factor {T/uqU) originates from the Bose occupation 
numbers of thermal excitations with energies of order ngt/, which are the most important in the perturbation theory. 
As shown below, the damping of excitations with energies Si, ^ tiqU is determined by both the Szepfalusy-Kondor 
(i' + 7 <-> 7') and Beliaev (z/ <-^- 7 + 7') processes, and can no longer be treated as Landau damping. 

The theory is generalized for the case of excitations in a spatially inhomogeneous (trapped) Bose-condensed gas. 
A new ingredient here is related to the inhomogeneous density profile of the condensate and the discrete structure of 
the excitation spectrum. We develop the method of calculating the self-energy functions and derive the equations for 
finding the wavefunctions and energies of the excitations (generalized Bogolyubov-De Gennes equations). 

The use of the theory is demonstrated by two examples. The first one concerns quasiclassical low-energy excitations 
of a trapped Bose-condensed gas in the Thomas- Fermi regime. The term "low-energy" assumes that the excitation 
energy is much smaller than the mean field interparticle interaction nQmU (nom is the maximum condensate density), 
and the quasiclassical character of the excitations requires the condition Si, hcu, where uj is the characteristic trap 
frequency. We consider anisotropic harmonic traps, where the discrete structure of the excitation spectrum is not 
important (see below and in ||2l[| ). On the contrary, the inhomogeneity of the condensate density profile has a crucial 
consequence for the damping rates and energy shifts of quasiclassical low-energy excitations. The most important 
turns out to be the boundary region of the condensate, where no{r)U ~ e^, |2^]. Therefore, the result for the damping 
rates and energy shifts is completely different from that in spatially homogeneous gases. 

Finally, we analyze the frequency shifts and damping of axially propagating sound waves in cylindrical Bose con- 
densates. As found, the nature of damping is similar to that in the case of phonons in spatially homogeneous Bose 
condensates. We show that the attenuation of axially propagating sound wave packets in the recent MIT experiment 
p3| can be well explained as a consequence of this damping. 



II. GENERAL EQUATIONS 



We consider a weakly interacting Bose-condensed gas confined in an external potential V{r). The grand canonical 
Hamiltonian of the gas can be written a,s H = Hq + Hi , where (hereinafter m = h = 1) 

Ho^ J dV#t(r) (^-^ + V{r) - /i^ *(r), (2) 

and the term 

Hi = j J d3r^t(r)*^(r)^'(r)#(r), (3) 

assumes a point interaction between atoms. The field operator of atoms 4'(r) can be represented as the sum of 
the above-condensate part ^P' and the condensate wavefunction 'So — (^) which is a c-number. As the interparticle 
interaction Hi contains both terms conserving the number of above-condensate particles and terms transferring two 



2 



above-condensate particles to the condensate (or two condensate particles to the above-condensate part) , the diagram 
technique should include both the normal Green function Q and the anomalous Green function (see, e.g. |p^ ). 




FIG. 1. The set of diagrams contributing to the normal self-energy S. Here a solid line with an arrow represents the normal 
Green function Q, solid line without an arrow corresponds to the anomalous Green function JT, white circle stands for the 
interaction vertex U and the black circle represents a sum of two white circles, one being a direct interaction and the other an 
exchange interaction. Dashed lines stand for the condensate wave function y/no. The self-energy part E"*" can be obtained by 
a time-reversal(i.e. the change t —> —t and p —> — p) of the graphs shown above. 




FIG. 2. The set of graphs contributing to the anomalous self-energy Ea. The notations are the same as for the Fig.l. 
The sums of the contributions of all irreducible diagrams will be represented by the normal (S) and anomalous 
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(Sq) self-energies (see Fig.l and Fig. 2). The former corresponds to the processes conserving the number of above- 
condensate particles, and the latter describes absorption (or emission) of two particles to (out of) the condensate. 
The Green function and self-energy operators satisfy Beliaev-Dyson equations [p^ , p^ 



g = G + G^g + G^a^, 
T = G+Y.+T + G+Y.ag, 



(4) 
(5) 



where the Green functions G and G^ describe forward and backward propagation of a particle characterized by the 
Hamiltonian Hq. 

We confine ourselves to the case of repulsive interaction between the atoms (a > 0). To develop the finite tem- 
perature perturbation theory for calculating dynamic properties and finding the excitation spectrum of a weakly 
interacting Bose-condensed gas we will use the non-equilibrium generalization |2j] of the Matsubara diagram tech- 
nique. In Eqs. (^),(|^) we perform an analytical continuation of the Matsubara frequencies Q — ^irTj (j is an integer 
number) to the upper half-plane, which corresponds to the replacement iQ e + iQ. Then, multiplying both sides of 
Eqs. (^) and (^) by G~^ and (G"*")"^, respectively, we arrive at the system of equations in the frequency-coordinate 
representation: 



eg{e;r,r') 
-eT{e; r, r') 



A 

T 

A 

T 



+ y(r)-/i + E(e) 
+ y(r)-/i + S+(£) 



g{e- r, r') -f ^a{e)T{e- r, r') + 5{y - r'), 
^(e;r, r') + ^a{e)g{e;v,v'). 



(6) 
(7) 



Here the action of the integral self-energy operators on the Green functions is written in a compact form 
J (Pr"Ti{e\ r, r")5(e; r", r') = S(e)C/(e; r, r') (and similar relations for the other combinations). Eqs. (^),(0) should be 
completed by a generalized Gross-Pitaevskii equation for the condensate wavefunction: 



-- + y(r)-Ai+(E 



^'o(r) = 



(8) 



and by the normalization condition 



dV(no(r) +n'(r)) = iV, 



where no(r) = |5'o(r)p is the condensate density, n'{r) = (^'^'^ {y)"^' [r)) is the density of above-condensate particles, 
and N the total number of particles in the gas. 

Eqs. (^,(0) can be solved by using the Bogolyubov transformation for the Green functions: 



g{e^] r, r') = u^(rX(r') + v^{y)vI{y'), 
T{e^]Y, r') = -u^{y)vI{y') - v^(r)u*(r'). 



where the index u stands for the set of quantum numbers, and the functions u^, Vi, satisfy generalized Bogolyubov-De 
Gennes equations 



£uUi,{y)-- 



A 

T 

A 

T 



+ y(r)-M+S(£,) 
+ y(r)-/i+E+(e,) 



-u^(r)-Sa(eJu,,(r), 
Vu{r)-T.a{e^)uAY) 



(9) 
(10) 



In the Bogolyubov-De Gennes approach only the terms bilinear in ^' operators are retained in the interaction 
Hamiltonian iJi, which assumes that the condensate density is much larger than the density of above-condensate 
particles. Then, the self-energy operators take the form 

E(£,r,r') = 2no[/<5(r-r'), (11) 
Ea(£,r,r') =noC/(5(r-r'). (12) 

The result of their action on the condensate wavefunction ^'o(r) and the functions u^y), Vu{y) is reduced to 
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S(£,)«'o(r) = / dV'S(£,,r,r')^'o(r') = 2no{r)U^o{r) 



and similar relations for the other combinations. Then, Eq.(@) becomes the ordinary Gross-Pitaevskii equation 



*o(r) = 0, 



and Eqs. (^),(0) are transformed to the ordinary Bogolyubov-De Gcnnes equations 



-j + V{v)-fi + 2no{r)U 



Uy(r) - no{r)Uv^{r), 



Taking into account Eq. ( ]l3| ) , in terms of the functions = Uy ±Vy these equations can be rewritten as 



(13) 

(14) 
(15) 

(16) 
(17) 



For a trapped Bose-condensed gas in the Thomas-Fermi regime, where /x ~ ?^omax 

U is much larger then the spacing 
between the trap levels, the kinetic energy term in Eq. (O) can be omitted and one has p5[ 



//i- V(r) 



(18) 



if the argument of the square root is positive and zero otherwise. For the low-energy excitations (e,y <C »omX/) 
of Thomas- Fermi condensates Eqs. (p^),(p7|) can be reduced to hydrodynamic equations obtained by Stringari |2^] 
and solved in the case of spherically symmetric harmonic potential V{r) and for some excitations in a cylindrically 
symmetric potential. An analytical method of solving Eqs. (p^),(p7[) (or the corresponding hydrodynamic equations) 
for the low-energy excitations of Thomas- Fermi condensates in an anisotropic harmonic potential V{y) has been 
developed in [|7|@. 

For a spatially homogeneous gas the generalized Gross-Pitaevskii equation (H) is equivalent to the Pines-Hugenholtz 
identity . In the Bogolyubov approach it simply gives /i = npC/, and Eqs. (p^),(|l5|) lead to the Bogolyubov spectrum 



(pV2)2 + noC/p2^ 



(19) 



where p is the momentum of the excitation. 

Under the condition rtg S> n', for which the Bogolyubov approach was originally developed, one can simply put no 
equal to the total density n in Eq.(n9|). For n' ^ uq, which can be the case a.t T ^ n, the dispersion law becomes 
essentially temperature dependent . In a spatially homogeneous gas the temperature dependence predominantly 
originates just from the presence of above condensate particles, with n' « n{T/Tcf^^ where = S.Sln^/^ is the 
BEG transition temperature. This leads to the replacement no ^ no + n' in Eq.(|Tl|) and gives /i = (no -I- 2n')U. 
The dispersion law will be still given by Eq.([l9|) in which the condensate density is now temperature dependent: 
no = n[l-(T/T,)3/2]. 



III. SPATIALLY HOMOGENEOUS BOSE-CONDENSED GAS 



In this section we present the results for the damping rates and energy shifts of elementary excitations in an infinitely 
large spatially homogeneous Bose-condensed gas. As one can see from Eqs.(p^,(|l5|), for finding the energy spectrum 
and wavefunctions of the excitations it is sufficient to calculate the self-energies S, and Sa. We will perform the 
calculations in the frequency-momentum representation and for physical transparency consider temperatures 

T > noU (20) 
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(the opposite hmiting case has been discussed by Popov ||T^ with regard to the phonon branch of the excitation 
spectrum). In the zero order approximation in the parameter (nQO^^^^ we have the weh-known mean field result: 
j](0) ^ 5^(0)+ ^ 2(no +n'(°))C/, = uqU, with = nfT/Tc)^/^ (see above). In this approach we obtain 

the Bogolyubov quasiparticle excitations with the spectrum dly), which we use in order to calculate the next order 

in [noa^Y^^. The latter is determined by the contribution of diagrams containing one quasiparticle loop Q (see 
Figures 1 and 2). Actually in this approach we represent the Hamiltonian as the sum of the (diagonalized) Bogolyubov 
Hamiltonian and the perturbation Hint originating from Hi (^ and containing the terms proportional to ^'o'l'''^ and 

Hint = U I dV[5'o(r)#'t(r){*'t(r) + *'(r)}*'(r) + (l/2)*'t(r)*'t(r)#'(r)4''(r)]. (21) 



Retaining only the temperature-dependent contributions, after laborious calculations for the normal self-energy we 
obtain E = S^") + S^^), where 

I](i)(P)=S"(F) + E'^(P), 

^"(^) - 2t/^no !4^An,^ n.) ^^^^^^ + ^^^^ " ^^^^^ + ^^'^^ 



(2^)3^ • e-Eq-eu 

2AqBk + BqBk — 4:BkCq + 2CqCk 



e + Ek+Eq 



8{Trnoa^)^/^T, (22) 



S'-(P) = 2U\ f,'lt'^ {2AqAk + 2AkBq + 2BkBq+4CkCq-AAkCq-4BqCk). (23) 

J [^7T ) E -rEq Ek 

Here P = {e, p}, k = q+p, Ep = p'^/2, Ep is given by Eq.(pj|), Uq is the equihbrium occupation number, Cp — naU /2Ep 
and Ap, Bp (±ep + Ep + noU)/2Ep. Similarly, the correction to the anomalous self-energy is given by 

v"^m 9r72 f d^Q t ^ f 2AkBq - 2AkCq - 2BqCk + 3CqCk 
J:^[P) = 2U no / j^-7^{nk +nq) ' 



(27r)3 \^ £-Sq-Ek 

2AqBk 2BkCq ~ 2AqCk + ?>CkCq\ _ 4(^„^^3)l/2r^ (24) 

E + Ek + Eq J 

j:iiP) =2U^no [-^-^hL^-^{2AkAq + 2BkBq + QCkCq~2ACq-2A,Ck-2BqCk-2BkCq). (25) 

J {ITT}-^ E + Eq~ Ek 

The resonant parts YJ" , SJ^ originate from the terms where one of the intermediate quasiparticles is created and 
another one annihilated, and the non-resonant parts S",!]" from the terms where both intermediate quasiparticles 
are created (annihilated). Temperature independent terms in the non-resonant parts, found by Beliaev |12[| , are 
omitted in Eqs.(|2|)-(|5l). 

Each of the self-energies (|2])-(p5h is singular at P — > and at least for small momenta the corrections become 
larger than the mean field values (p|). Nevertheless, keeping in mind that any physical quantity is determined by the 
combinations of the self-energies, which do not contain the infra-red singularities, we will still treat E*^^^ and si^-* as 
perturbations. 

For a spatially homogeneous gas the Pines-Hugenholtz identity /i — (S(P) — T,a{P))\p^Q gives the first order 
correction to the chemical potential 

= [i = 12(^3)1/2^^ (26) 

and the relation between no and the chemical potential, /i = noU — /3^Vno, coincides with that found by Popov 
[Q. The u,v functions in generalized Bogolyubov-De Gennes equations (g), ( p^ ) can be written as Upexp(ipr) and 
?;pexp(ipr), and in terms of the functions — Up±Vp these equations take the form 



(e-5-(P))/-= (^^ + ^1(P))/;, 



(27) 



S-iP))f^ ={4 + 2noU + S+iP) ] /-, (28) 
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where 



S^^ ^ +/3V^±E„ (29) 

S- = . (30) 



Considering the terms S , S±. in Eqs. (|27|) , (|28D as small perturbations we put e = Sp in the expressions for these 

quantities, following from Eqs. (p2[)-(p5|). Then, solving Eqs. (|27|) , (p8|) , for the excitation energy we obtain e = ep+s^\ 
where 

(31) 



p-(ep,p) 



As S*-^-* and si"'^'' are complex, the correction to the excitation energy has both a real and an imaginary part: 
Ep^'' = Sep — iVp. The former gives the energy shift, and the latter is responsible for damping of the excitations. Under 



the condition T ^ not/ a straightforward calculation of Eq.(plD on the basis of Eqs. (22)-(^5|) and (|29|),(pOD, for the 
phonon branch of the spectrum (sp <ti jiqU) yields 

fep«-7£p^(noa3)i/2, (32) 
noU 

Tp=ep'-^{n,a'f'\ (33) 

It is important to emphasize that in this case the energy shift is determined by both resonant and non-resonant terms, 
whereas the damping rate Tp is described solely by the resonant contributions. This type of damping originates from 
quasi-resonant scattering of thermal excitations from a given excitation (Landau damping) and is absent at T = 0. 
Both the energy shift and the damping rate are determined by the interaction of a given excitation with intermediate 
quasiparticles having energies Sq ^ uqU . The damping rate Tp (^ ) coincides with that found in recent contributions 



20 -22 1 and contains a slight numerical difference from the earlier Szepfalusy-Kondor result |Q. The energy shift 
for the phonon branch of the spectrum was also calculated in . In the latter work the expansion of the self-energy 
functions near the point e = Ep was used and formally divergent integrals were canceling each other in the final 
expression for the energy shift, which have led to the result by approximately factor 6 smaller then the shift ( |3^ ) 
obtained by the exact integration. 

Eqs.(p2|),(|3^) clearly show that the small parameter of the theory is (T/no?7)(fT.oa'^)^^^ ^ 1 (see Eq.(|l|)), whereas in 
the zero temperature approach |^ the small parameter is (nga^)^/^ <C 1. The presence of an additional large factor 
T/uqU at finite temperatures T 3> uqU originates from the Bose enhancement diagrams containing one quasiparticle 
loop: Compared to the zero-temperature case the contribution of each of these diagrams is multiplied by the Bose 
factor Uq — [exp(eg/T) — 1]~^ (or l + Uq). As the most important is the contribution of intermediate quasiparticles 
with energies Eg ~ noU, for T 3> riQU the Bose factor riq ^ T/uqU. The criterion similar to Eq.(0) was found by Popov 
P,p^ as the condition which allows one to use the mean field approach at finite temperatures and to renormalize the 
theory for reaching beyond this approach. 

Remarkably, the criterion (|l|) is fulfilled even at temperatures very close to Tc- For AT ^ Tc — T Tc we have 
no ~ n/ST /Tc, and Eq.(|l|) gives AT /T^ ^ (noa^)^/'^, which coincides with the well known Ginzburg criterion ||2^ for 
the absence of critical fluctuations. The criterion also ensures that the main contribution to the damping rate 
originates from the interaction of a given excitation with thermal excitations through the condensate, i.e., from the 



first term in Hi-at (21) 



At energies Ep ^ noli (but still Ep <C T) the results for the energy shift and damping rate, following from Eq.(^, 
can be represented in the form 

5ep = EpFJ^)^{noa')y\ (34) 
\noLJ / noU 

Tp^EpFJ^)^{ma^r/\ (35) 
\noU J noU 

The functions Fi and F2 have been calculated numerically. For Ep <ti uqU we have Fi k. —7, F2 = 37r'^/^/4 and 
Eqs. (^), (^H) coincide with Eqs. (^), (^). The dependence of 5Ep and Fp on the excitation energy is presented in 
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Fig. 3. For Sp ^ uqU the damping rate increases with Ep and reaches its maximum at Sp ^ IQuqU. Further increase of Sp 
leads to slowly decreasing Pp. For single-particle excitations {Sp = tiqU) we obtain F2{ep/noU) ~ {noU /Sp)^^'^ . 

Accordingly, the damping rate can be written as Tp {T/ep)noaVp, where a — Stto^ is the elastic cross section, and 
Vp the particle velocity. This damping rate exceeds the Beliaev temperature-independent term naVp even at T close 
to Tc, if AT = Tc — T ^ Ep. In contrast to the phonon branch of the spectrum, for Sp ^ uqU the damping is provided 
by both the Szepfalusy-Kondor (v 7') and Beliaev (z^ <-> 7 -I- 7') processes and, hence, can no longer be treated 

as Landau damping. The small parameter of the theory is still given by Eq.(|^), since even at Sp ^ noU the energy 
of at least one of the thermal excitations is of order uqU. 




X 



FIG. 3. The damping rate Tp (solid line) and the energy shift Sep (dashed line) versus {sp/noU). Both Tp and Ssp are given 
in the units of r(noa^)^''^. 



The energy shift for Sp ^ noU is negative. The modulus of the shift increases with Sp and reaches its maximum 
at £p « AnoU. The further increase of Sp decreases |fep|. The latter is equal to zero for Sp w GOnoU, and becomes 
positive at larger Sp. 

The above results for the damping rate and energy shift of a given excitation are obtained in the so called coUisionless 
regime: We assume that the De Broglie wavelength of the excitation, 1/p, is much larger than the mean free path of 
the thermal quasiparticles with energies ~ uqU , which are mostly responsible for the damping and shifts. It is also 
assumed that the excitation energy Sp greatly exceeds the damping rate of these thermal excitations. The latter is 

of order T{noa^)^^'^ (see Fig. 3), and for £p > noU the two requirements of the coUisionless regime are well satisfied 
under condition (|l|). In the phonon branch of the excitation spectrum {sp <C ng?/) these requirements are equivalent 
to each other, and the coUisionless criterion can be simply written as 

Sp > T{naay/\ (36) 

As clearly seen, in the phonon branch one can always find excitations which do not satisfy Eq.(|3^) and, hence, require 
a hydrodynamic description with regard to their damping rates and energy shifts. 

The coUisionless criterion (p6 ^provides an additional argument on support of the above used perturbative ap- 
proach for solving Eqs. (p^), (2m. Under condition (|36|) the term 5" < T{noa^y/^ < Cp, the term < 
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ep(T/noC/)(noa3)i/2 < g^^^ ^^le term 5*+ < {uqU /ep)T{noa^)^/'^ < UqU. 

The non- mean- field shift is actually the shift of the excitation energy Sp at a given condensate density uq. On 
the other hand, Sp is determined by the Bogolyubov dispersion law (p^, with the temperature-dependent condensate 
density rio(T), and, hence, is temperature-dependent by itself. Therefore, at a given T one will also have the mean- 
field temperature-dependent energy shift Se™^ = £p{T) — ep(0). As the condensate density decreases with increasing 

temperature, 5e™^ is always negative. For T 3> noU it greatly exceeds the above calculated shift 5ep at any p. The 
ratio (fc™^/(5ep) decreases with temperature, but even for n' <C tt-o one has 



5ef = -ep{{)) ^^^^^1^, (37) 



and {5ef/5ep) - (T/noC/)!/^ > 1. 



IV. SPATIALLY INHOMOGENEOUS BOSE-CONDENSED GAS 

We now generalize the above obtained results for the case of elementary excitations in a spatially inhomogeneous 
(trapped) Bose-condensed gas. As already mentioned in the introduction, a new ingredient here is related to the 
inhomogeneous density profile of the condensate and the discrete structure of the excitation spectrum. This requires 
us to develop a new method of calculating the self-energy functions in generalized Bogolyubov-De Gennes equations 
(^, (p^. The self-energy operators in these equations are the sums of the zeroth and first order terms: 

E,(£, r, r') - 7io(r)c7<5(r ~ r') + (38) 
I](e, r, r') = 2(no(r) + n"^''^)U6{r - r') + S^^), (39) 

and a similar relation for S+. At temperatures T 3> nQ„iU the zero order value of the above-condensate density in 
the condensate spatial region is coordinate independent and equal to the above-condensate density in the ideal gas 
approach: n'(°)(T) = 2.6{T/2nf/^. On the contrary, the self-energies % Ea depend explicitly on the condensate 
density. Due to the discrete structure of the energy spectrum of excitations the expressions for these quantities should 
be written in the form of sums over the discrete states of intermediate quasiparticles 7,7'. In the frequency-coordinate 
representation we have 

S(£,, r, r') = E"(£,, r, r') + £'■(£,, r, r'), 

I]"(e., r, r') ^ 2noU' Y.^n, + ny) r 2%W"7(i-Ov(r)v(rO + u,(r)^.,(rO.,,(r).,,(rO , 
—4:U-y{r)uy{r')uy'{r)vy (r) + 2u-y(r)w^(r')M^/(r)w^/(r') 



e + e-y + Sj' 
2u^{r)v^{r')uy{r)vy{r') \ _ ^^3.1/2 



s + e-y + e-y' 



E'^(e„r,r') = 2nof>'^ 



-8{noTTa^y'^TS{r-r'), (40) 



7,7' 



{2u~^{v)u^{y')u^i {y)u^i {v') + 2u^{r)u^{v')v^i{v)v^i{r') + 2v^{v)v~^{y')v^i {r)v^: {v') + 
4u^(r)t;..y(r')-u^'(r)w^'(r') — 4M^(r)u^(r')u^/ (r)u^/ (r') — '^u^{r)v^{v')v^i {r)v^i {v')) . (41) 

S„ = E^(e„r,r') + S',(£.,r,r'), 
E;'(£^, r, r') = 2nQU'^Y^{n^ + ) 
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j:iie,,r,r')^2noU^J2 



S I £ -y I S^'-y' 



(42) 



1 



7:7' 



-7 



(2u^(r)u-y(r')u,y/(r)u,y' (r') + 2vj{r)vj{r')vj/{r)vji{r') + 6uj{r)vj{r')uji {r)vji {r') — 
2u^{r)u^{r')uj' {r)v^/ {r') — 2u^{r)v^{r')u^/{r)u^/{r') — 2u^{r)v^{r')vY {r)v^/ {r') — 
2v-y{r)v^{r')uj/{r)v^' {r')). 



(43) 



As we saw in the previous section, in the spatially homogeneous case all physical quantities are determined by the 
contribution to the self-energy functions E'^^\ eI^-* from intermediate quasiparticles with energies of order the mean 
field interaction between particles. The same holds for a spatially inhomogeneous (trapped) Bose-condensed gas in 
the Thomas-Fermi regime, where the mean field interaction nQmU greatly exceeds the level spacing in the trapping 
potential. The intermediate quasiparticles with energies of order nomU are essentially quasiclassical. With regard 
to the integral operator (S^^^ — I]i^^)e_^o in the generalized Gross-Pitaevskii equation (||), which is solely determined 
by non-resonant contributions, this immediately allows one to replace the summation over the discrete intermediate 
states by integration. The kernel of this integral operator varies at distances |r — r'| of order the correlation length 



^cor — ^/ynoU which is much smaller than the characteristic size of the condensate. Therefore, the result of the 
operator action on the condensate wavefunction can be written in the local density approximation and, hence, should 
rely on Eq.(pq) with coordinate-dependent condensate density no(r): 



6^0 



eW-EW] *o(r)= /dV[E«(£,r,r')-E«(£,r,r')l ^oir') ^ -f3^lir). 



(44) 



This result can be easily obtained from Eqs. (^), (^), where one should put = 0, neglect the difference between 
e-y and e-yi, and make a summation over 7'. Replacing the summation over 7 by integration one should also take into 
account that for quasiclassical excitations the functions can be represented in the form 



/£2 + (no(r)f/)2-no(r)[/ 



Tl/2 



AW, 



(45) 



where |/.y(r)p is the ratio of the local to total density of states for Bogolyubov quasiparticles of a given symmetry, 
described by the classical Hamiltonian 



if(p,r) = ^(pV2)2+no(r)!7p2. 
On the basis of Eq. (ji^ ) we obtain the generalized Gross-Pitaevskii equation in the form 

-j + V{r) -fi + U\^o\^ - /3*o) *o = 0, 



(46) 



(47) 



where p, = 11 — 2n'^'^\T)U is coordinate independent. Compared to the ordinary Gross-Pitaevskii equation (p^), 
Eq.(^7|) contains an extra term [2n'^^^U — /3^'o]5'o hi the Ihs. One can easily check that Eq.(^) coincides with the 
equation 
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obtained by averaging the non-linear Schrodinger equation for the field operator. As mentioned above, for T ^ uqU 
the above-condensate density n' — (^P'tvl/') in the condensate spatial region is mainly determined by the coordinate- 
independent (ideal gas) contribution n''^°) (T) . The coordinate-dependent correction to the above-condensate density, 
n'(^)(r) = (5''t(r)^''(r)) — turns out to be equal to the anomalous average (5'(r)5'(r)): 



n'^^\v) = (*(r)1'(r)) 



3J7 



Accordingly, the quantity [2n'^^'^U - (3^g\^q ^ [/(2(*'t*') + (*'*')). 

Taking advantage of Eqs (|3^), ( p9[ ) and (p7|), the generalized Bogolyubov-De Gennes equations (||), ( p^ are reduced 

to 



- 5-)/-(r) = ( + ||£ + /+(r) 
(s. ~ 5-)/+(r) = ( -| + 1^ + 2|*opt^ + S] 



(48) 
(49) 



are given by Eqs. (29), (pOl), and is the exact value of the excitation energy. Eqs. M 



where the quantities S^, S~ 

( ^ ) and ( ^ ) represent a complete set of equations for finding the energy shifts and damping rates of the elementary 
excitations. 

A precise calculation of the self-energy functions in Eqs. (^), ( ^ ) depends on the value of e^, and on the trapping 
geometry. In this section we will make general statements on how the calculation can be performed. In most of the 
cases (except the case of the lowest excitations with zero orbital angular momentum in spherically symmetric traps) 
the characteristic time scale in the self-energy operators, l/ffzy, is much smaller than the inverse level spacing in the 
trap. Therefore, the summation over the discrete intermediate states can be replaced by integration. This is a direct 
consequence of the general statement that the time-dependent discrete Fourier sum can be replaced by its integral 
representation at times much smaller than the inverse frequency spacing (see e.g. , p^). 

The kernels of the non-resonant parts of the self-energy operators, (Eof) and (p2), vary at distances |r — r'| which do 



not exceed the correlation length Zcor = l/y noU. As Zcor is much smaller than the characteristic size of the condensate, 
the non-resonant parts of the self-energies can be calculated in the local density approximation. One can easily find 
from Eqs. (E^), (^) that they lead to the result of Eqs. (p2|), (24) (with hq replaced by the coordinate-dependent 
density no{rJ), multiplied by S{r — r'). For the lowest excitations one should put k = q in Eqs. (p^) , ([2^) , and for 
quasiclassical excitations take p from the Bogolyubov dispersion law H{p, r) = e^. 

The calculation of the resonant contributions to the self-energies is more subtle. Using Eq.(^) for the functions 
Ui,,Vi, — in Eqs. (^l]), (43), one can see that all resonant contributions contain the quantity 



Q(r,r') = E 

7' 



/^(r)/^(r')/y(r)P(r') 



iO 



Writing {si, + 



iO) ^ as the integral over time exp{i(e,y + — e^i -I- iO)i}, we obtain 

r-oo 

Q{r,r') = i / dteyip {iSi,t)K^{r,r' ,t), 
Jq 



where the quantum-mechanical correlation function 



K^{r, r', t)^J2 /7(i-)/7(r')/7'(r)/r (r') exp {i{e^ - + iO)t} . 



(50) 



(51) 



We will turn from the integration over the quantum states 7' of the quasiclassical thermal excitations to the integration 
along the classical trajectories of motion of Bogolyubov-type quasiparticles in the trap. Following a general method 
(see [^S]-^), employed in for the damping of low-energy excitations, we obtain 



K,{r, r',t)^ g~' J 6{r" - r)d{r';{t) - r')6{e^ - H{p, r")) 



87r3 ' 



(52) 



where Yp(t) is the coordinate of the classical trajectory with initial momentum p and coordinate r. Eq.(|5^) will be 
used in the next sections where we demonstrate the facilities of the theory. 
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Concluding this section, we emphasize the key role of harmonicity of the trapping potential for temperature- 
dependent energy shifts of the excitations. As mentioned in the previous section, in the spatially homogeneous case 
at a given temperature the non-mean-field shift is much smaller than the shift Se"^^ appearing in the mean field 
approach simply due to the temperature dependence of the condensate density in the Bogolyubov dispersion law 
(p^). For the Thomas-Fermi condensate in a harmonic confining potential the situation is different. In this case the 
spectrum of low-energy (e^ <ti nomU) excitations is independent of the mean field interparticle interaction nomU 



(chemical potential) and the condensate density profile 



8| . Hence, the temperature-dependent energy shifts can 



only appear due to non-Thomas-Fermi corrections. For finding these corrections one should use the mean field self- 
energies I]a(£, r, r') — no(r)[/(5(r — r'), r, r') = 2{no{r) + n'^^") )U S{r — r') , where the only difference from the T ~ 
case is related to the presence of above-condensate particles in the condensate spatial region at finite T through the 
coordinate-independent term 2n'(*'^[/ in E. Then Eqs. (|4^), ( p9[ ) take the form of ordinary Bogolyubov-De Gennes 
equations ([l6|) , (|l7|) , and Eq.(^7|) becomes the ordinary Gross-Pitaevskii equation (|3|), with the chemical potential /i 
replaced by fl. The latter circumstance changes the condensate wavefunction compared to that at T = and ensures 
the temperature dependence of ^'q. Accordingly, the excitation energies Si, in Eqs. (p^,([l7|) also become temperature 
dependent. This type of approach, which for a spatially homogeneous gas would immediately lead to the result of Lee 
and Yang ||^ , has been used in recent numerical calculations of the energy shifts of the lowest quadrupole excitations 
in spherically symmetric |^o| and cylindrically symmetric |^,^ harmonic traps. The presence of the coordinate- 
dependent part of the above-condensate density, n'^^-'(r), in these calculations is not adequate, since the anomalous 
average equal to this part was omitted and equations for the excitations did not contain the corrections to the self- 
energies, also proportional to {noa^Y^^. However, at T ^ noU, where rt'^^^ ^ the coordinate-dependent part 
n''^^-'(r) as itself should not significantly influence the result, and the calculations should actually demonstrate 
how important are the mean field non-Thomas-Fermi effects. The results of [ |3l] , |32| | show the absence of energy shifts 
of the excitations at temperatures T < 0.6Tc in the JILA experiment [Q and in this sense agree with the experimental 
data, but do not describe the upward and downward shifts of the excitation energies, observed experimentally at higher 
temperatures (in this respect it is worth mentioning that the calculations |]33|] performed for the thermal cloud in the 
hydrodynamic regime agree surprisingly well with the experiment Q). On the other hand, the calculation shows a 
downward shift of the energy of the lowest quadrupole excitation with increasing temperature in the conditions of the 
MIT experiment This is consistent with the experimental data and indicates that for not very small Thomas-Fermi 
parameter Lu/nQmU the mean field non-Thomas- Fermi effects can be important for temperature-dependent shifts of 
the lowest excitations. ^ 

Below we will assume a sufficiently small Thomas- Fermi parameter uj/nQ„iU and demonstrate the use of the theory 
by the examples where the influence of non- Thomas-Fermi effects on the energy shifts of the excitations is not 
important. 



V. QUASICLASSICAL EXCITATIONS IN A TRAPPED BOSE-CONDENSED GAS 

We will discuss the Thomas-Fermi condensates in a harmonic confining potential on the basis of Eqs. (^7|), (|4^), 
(^). Neglecting the kinetic energy term in Eq.(^), we arrive at a quadratic equation for "^Iq. Expanding the solution 
of this equation in powers of f3 and retaining only the terms independent of f3 and the terms linear in /3, for the 
condensate density we obtain 

noir)^fio{r)+Py/f^/U, (53) 

where no(r) — {jl — V{y))/U is the density of the Thomas-Fermi condensate in the ordinary mean field approach. 

We first consider the damping and energy shifts of quasiclassical (e^ ^ to) low-energy excitations of a trapped 
Thomas- Fermi condensate, i.e., the quasiclassical excitations with energies much smaller than the mean field interac- 
tion between particles nomU. In this case the terms in Eqs. (^) and (^9|), originating from the kinetic energy of the 
condensate, can be omitted from the very beginning ^t\ . Then, using Eq.(|5^) and treating the terms containing S~ 
and as perturbations, we obtain e^, = + e'i}\ where el?-' is the excitation energy in the mean field approach, 
and the correction to the excitation energy = Ssi, — iVi, is given by the relation 

ei'^ = {f-\S-\f+) + \ + 2p^\f-) + . (54) 

Here /^^°^ are the zero-order wavefunctions of the excitations, determined by the ordinary Bogolyubov-De Gennes 
equations ([T6|), (|l7|), with Sy — Si/ 
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In the case of quasiclassical excitations also the kernels of resonant parts of integral operators in Eq. ( p4| ) vary on 
a distance scale |r — r'| which does not exceed the correlation length /cor- This can be already seen from Eqs. (|5(j), 
(p^): The characteristic time scale l/sv in Eq.(^0|) is much shorter than u)~^ and important is only a small part of the 
classical trajectory, where the condensate density is practically constant and Vp{t) = r + vi, with v — dH/dp. The 
correlation length /cor is not only much smaller than the size of the condensate, but also smaller than the width of the 
boundary region of the condensate, where noU ^ e„. Therefore, the action of all integral operators on the functions 
/^(°) in Eq.® can be calculated in the local density approximation. Accordingly, for each of these operators one 
can use the quantity following from Eqs. (|2^)-(|2^), with no = no(r) and p from the Bogolyubov dispersion law 
H{p, r) = e^. Then, using Eqs. (^) we can express the energy shift Se^, and the damping rate through the energy 
shift fc,^/j(r) and damping rate T^}i(r) of the excitation of energy e^, in a spatially homogeneous Bose-condensed gas 
with the condensate density equal to no(r): 



Se^ldMMr)n5e.H{r)+^^d^Sl }, (55) 



= / dV|/,(r)|2r,^(r). (56) 



The second term in the integrand of Eq. ( |55|) originates from the temperature dependence of the shape of the condensate 
wavefunction. For any ratio s^/noir)!! this positive term dominates over the negative term &[^?i(r). The latter 
circumstance can be easily established from the results for Se,^h in Fig. 3. Thus, for quasiclassical low-energy excitations 
the energy shift Sei, will be always positive, irrespective of the trapping geometry and the symmetry of the excitation. 

We confine ourselves to the case of cylindrical symmetry, where for the states with zero angular momentum one 
finds 

l/.(r)P = ^ , (57) 



nlph log {2jl/Eu)p\J eI + (no(r)C/)2 



with Ip = (2/i/[x)p)^/^, Iz = {2fi/ijJzY/'^ being the characteristic size of the condensate in the radial and axial direction, 
cjp, ljJz the radial and axial frequencies, and p the radial coordinate. The main contribution to the integral in Eq.(|55|) 
comes from the boundary region of the condensate, where no(r) ~ e^- From Eqs. (p8[),(p9|) one can easily see that in 
this region the possibility to omit the non-Thomas- Fermi effects originating from the kinetic energy of the condensate 
requires the condition Ei, '3> w'^^^p}/^. This condition ensures that the characteristic width of the boundary region 
greatly exceeds the excitation wavelength, and we arrive at the following relations for the energy shifts and damping 
rates of the excitations: 

fe. « ,J,, A non.a^f'\ (58) 

V M log 

,J,, A n,ma'Y''. (59) 
p log {2p/Ey) 

It is important to emphasize that in the boundary region of the condensate, responsible for the energy shifts 
and damping rates of the quasiclassical excitations, the quantities S~^^\ S'^^^\ and si^^ are determined by the 
contribution of intermediate quasiparticles which have energies comparable with E,y. Moreover, in this spatial region 
the quasiparticle energies are of order the local mean field interparticle interaction. As a consequence, the energy shift 
(5£j/ ( p8| ) and the damping rate Ty ( |59| ) are practically independent of the condensate density profile. For the same 
reason the damping rate is determined by both the Szepfalusy-Kondor and Beliaev damping processes. Therefore, 
similarly to the damping of excitations with energies Ep ^ noU in a spatially homogeneous gas, the damping of 
quasiclassical low-energy excitations of a trapped Bose-condensed gas can no longer be treated as Landau damping. 



VI. SOUND WAVES IN CYLINDRICAL BOSE CONDENSATES 



The derivation of Eqs. (^8]), (59) assumes that the motion of the excitation v is quasiclassical for all degrees of 
freedom. We now turn to the condensate excitations in cigar-shaped cylindrical traps, which are quasiclassical only in 
the axial direction and correspond to the lowest modes of the radial motion. We will consider low-energy excitations 
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(ffy <C nomU), i.e., the excitations with the axial wavelength much larger than the correlation length ?cor- In the 
recent MIT experiment ||2^ localized excitations of this type were created in the center of the trap by modifying the 
trapping potential using the dipole force of a focused off-resonant laser beam. Then, a wave packet traveling along 
the axis of the cylindrical trap (axially propagating sound wave) was observed. In the mean field approach the sound 
waves propagating in an infinitely long (axially homogeneous) cylindrical Bose condensate have been discussed in 

For revealing the key features of the non-mean-field effects (damping and the change of the sound velocity) we 
confine ourselves to the same trapping geometry. With regard to realistic cylindrical traps this will be a good 
approach if the mean free path of sound waves is smaller than the characteristic axial size of the sample. As found 
in pst , for axially propagating sound waves radial oscillations of the condensate are absent, and the wavefunctions 

/± = i2Mp)U/ek)^^/^fk, with 

fk = exp (ikz) (60) 

and k being the axial momentum. The dispersion law 

Ek = ck (61) 

is characterized by the sound velocity equal to (riomf//2)^/^, where nom — fi/U is the maximum density of the 
Thomas-Fermi condensate in the ordinary mean field approach. 

It should be noted from the very beginning that, according to Eq.(|5^), fiom is related to the corrected value of the 
maximum condensate density no™ as no™ = no™ — (/3/J7)^nom- Therefore, being interested in the sound velocity at 



a given value of the maximum condensate density, one should substitute this expression to Eq.(61). This immediately 
changes the sound velocity to 

c= (no™C//2)i/2 (62) 
in the leading term (^Tj) of the dispersion law and provides a contribution to the frequency shift of the sound wave 

6T 

5ek = -ek ^(7rno™a3)i/2. (63) 

The damping rate and other contributions to the frequency shift can be found directly from Eq.(|5^) by using 
the wavefunctions fk (|6^). The intermediate quasiparticles giving the main contribution to the damping rate and 
frequency shift have energies e nomU, i.e. much larger than the frequency of the considered sound wave, Sk (see 
below). Therefore, similarly to the case of phonons in a spatially homogeneous condensate, the non-resonant terms 
(^), (^ ) contribute only to the frequency shift. As already mentioned above, the characteristic distance scale \p — p'\ 
in the kernels of the self-energies (^0|), (^) is of order the correlation length Zcor, and the sum of their contributions to 
the frequency shift, (5e^?, can be calculated by using the local density approximation for the action of the self-energy 
operators on the functions f^. As a result, we express &^ through the non-resonant part of the energy shift S^f^ip) 
in the spatially homogeneous condensate with the condensate density no(p): 



For Ek ^ noU one can directly find from Eqs. (p2|), ( p4| ) that S^,^ = £k{T /noU){nnoa'^)^^'^ . Then we obtain 



fc^-£fc ^(7rno™a3)i/2. (64) 

nomU 

The resonant terms (|l]), ( ^ ) contribute to both the frequency shift and damping rate. This means that the latter 
is determined by the Szepfalusy-Kondor scattering processes and, since the characteristic energies of intermediate 
quasiparticles are much larger than Ek, can be treated as Landau damping. The resonant contributions to the 
frequency shift and damping rate can not be found in the local density approximation, as the characteristic distance 
scale \p — p'\ in the kernels of the self-energy operators in Eq.(p^ is of order the radial size of the condensate. For 
finding these contributions one has to substitute the resonant parts of the self-energies, (E^-(E3), to Eq.(|53) and. 
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by using Eqs.(p0[)-(p2[), turn from summation over quasiclassical states 7, 7' of intermediate quasiparticles to the 
integration along classical trajectories of their motion. Then, a direct calculation yields (cf. [^ ) 



e^U f dn f°° f 

J d^^-^ dtexpiiskt) J $fc^(r)$^^(r(p,i))5(£^-i7(p,r)) 



(27r)3 ' 



(65) 



where r(p, t) is the classical trajectory starting at the phase space points (r, p) on the (hyper)surface of constant 
energy Sj, <&fc^(r) = fk{z)F^{p), and 



2e. 



+ {m{p)Uf-m{p)U^e-'^ + {m{p)U? 



Generally speaking, the integration in Eq.(|65|) is a tedious task as it requires a full knowledge of the classical 
trajectories on a time scale ~ l/cfc. This is also the case in the idealized cylindrical trap, because of coupling 
between the radial and axial degrees of freedom. We will rely on the approach which assumes a fast radial motion of 
quasiparticles compared to their motion in the axial direction and, hence, requires the frequency of the sound wave, 
Ek, significantly smaller than the radial frequency ojp. Then on a time scale ~ l/e^ the quasiparticles with energies 
^ tiqU (which are the most important for the energy shifts and damping of the sound wave) oscillate many times in 
the radial direction, whereas their axial variables z{p, t), Pz{p, t) only slightly change and, hence, can be adiabatically 
separated from the fast radial variables p(p, i), Pp(p, i). In this case it is convenient to integrate Eq,(p5|) over dpp 
and, using Eq. (|60|) , represent it in the form 

51 = i^f2 j ^f^'^^'r^ J dtexpiiskt) J pdpdp.dz^^^MM^k^ exp {i{z - z{p,t))}, (66) 



+ (n„(p)[/)2 



where the integration is performed over the entire classically accessible region of the phase space. 

Since z(p,t) is close to z, in the exponent of the integrand we can write z(p,t) — z = v^, where the axial velocity 
Vz is obtained from the exact Hamiltonian equations of motion by averaging over the fast radial variables: — 
{dH{p,r)/dpz)p- For the classical radial motion (e-y ^ uip) the averaging procedure simply reduces to integration 
over dp under the condition iJ(p, r) = e-y at fixed values of e^, Pz and z, with the weight proportional to the local 
density of states for the radial motion: 



where g = /(e£ + (no(p)?7)^) ^^^2npdp. Finally, averaging the function F^{p{t)) over the fast radial variables and 



integrating over dt in Eq.(66), we obtain 



£ 
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U f , e^dn^ f , , {F{p)); 



61 = -1^ / de^^^—^ I dzdpzg ^' (67) 

SnHlJ ^ dsj J " Ek-PzVz+iQ 



The resonant contribution to the frequency shift, given by the real part of Eq.(|67|), after the integration proves to be 
(5e^ ~ —2.3ek{T/nQrnU){nomO'^y^'^- The sum of this quantity with the non-resonant term (|6^) and Ssk ( |63| ) leads to 
the frequency shift of the sound wave 

fefc « -4£fe(no™a3)i/'^^. (68) 
nomU 

The imaginary part of Eq. ( |67| ) gives the damping rate 

Ffe = 8.6efe(no™a3)V2^1^. (69) 

nOrnU 



Except for the numerical coefficients, Eqs,(|68|), (|69|) are similar to Eq.(p2D, ( |33D for the damping rate and energy 
shift of phonons in a spatially homogeneous Bose condensate. This a consequence of the fact that the condensate 
boundary region practically does not contribute to the damping rate and frequency shift of axially propagating sound 
waves, in contrast to the case of excitations quasiclassical for both axial and radial degrees of freedom. 
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In the MIT experiment [g3| the characteristic spatial size of created locahzed excitations was A « 20 /^m and, 
accordingly, so was the initial wavelength of propagating sound. According to the experimental data, the propagating 
pulse died out during 25 ms, and after that only the lowest quadrupole excitation characterized by a much loner 
damping time (~ 300 ms) was observed. We believe that the attenuation of axially propagating sound in the MIT 
experiment p3t on the time scale of 25 ms can be well explained as a consequence of damping. The characteristic 

frequency of the waves in the packet can be estimated as Si, ~ '^t^\J nomU /2/A. Then Eq.(p9|) gives the damping rate 
independent of the condensate density fiom- In the MIT experiment the temperature T «0.5 /iK was roughly only 
twice as large as nomU , which decreases the damping rate by approximately 20% compared to that given by Eq.(|69|). 
In these conditions we obtain a characteristic damping time of 15 ms, relatively close to the measured value. 

The relative change of the sound velocity, 5c/c = Ssk/sk, increases with decreasing condensate density n^m- How- 
ever, even at the lowest densities of the MIT experiment (?T-om ~ 10^^ cm~^) the quantity 5c/ c does not exceed 
~ 5% and is practically invisible. 
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